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We work out the magnetization and susceptibility of Heisenberg- and XXZ-model 
antiferromagnet spin- 1/2 systems in D dimensions under a rigorous constraint of 



single particle site occupancy. Quantum fluctuations are taken into account up to 
the first order in a loop expansion beyond the Neel state mean field solution. We 



discuss the results, their validity in the vicinity of the critical point and compare 



them with the results obtained by means of a spin wave approach. 
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Recent work on quantum spin systems discuss the possible existence of spin liquid states 
and in two-dimensional space dimensions the competition or phase transition between spin 



liquid states and an anti 



! Heisenberg type systems 

S3 
O 
O 



erromagnetic Neel state which is naturally expected to describe 



HQS 



show an antiferromagnetic phase 



. It is also known that undoped superconducting systems 



In the following we focus our attention on a Neel phase description of quantum spin 
systems described by Heisenberg models. More precisely we present below a detailed study 
of the magnetization and the parallel magnetic susceptibility of Heisenberg antiferromagnetic 
spin-1/2 systems on D-dimensional lattices at finite temperature. The aim of the work is 
the study of the physical pertinence of the Neel state ansatz as a mean-field approximation 
in the temperature interval < T < T c where T c is the critical temperature. In order to get 
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a precise answer to this point we work out the quantum fluctuation contributions beyond 
the mean-field approximation under the constraint of strict single site occupancy which 
allows to avoid a Lagrange multiplier approximation [10]. The results are also extended to 
anisotropic XXZ systems and compared to those obtained in the framework of the spin 
wave approach. 

The paper is organized as follows. In section II we present the derivation of the partition 
function under the single particle site occupation constraint. The mean-field and first order 
loop expansion term contributions are derived in section III. In section IV we determine 
the magnetization and the magnetic susceptibility, and discuss the results obtained at the 
different levels of approximation. Comments are presented and conclusions are drawn in 
section V. Details of calculations are presented in the appendix, section VI. 



II. FERMIONIZATION OF THE HEISENBERG MODEL AND THE 

PARTITION FUNCTION 

The Heisenberg antiferromagnet Hamiltonian (HAFM) in the presence of a local magnetic 
field Bi reads 

H = ~\ E J M + E*-3 w 

<i,j> i 

where Jjj < and the sums in the first term run over nearest-neighbour sites < i,j > on a 
D-dimensional hypercubic lattice. 

The S = 1/2 spin vector operators are expressed in terms of fermionic creation and 
annihilation operators {fi\,f} a } 



Si — fiaPa\fi\ (2) 

where the a a \ vector components are Pauli matrices. 

The transformation is rigorous if Y2 a flafia — 1- The Fock space constructed with the 
fermionic operators /, f> is not in bijective correspondence with the Hilbert space of the 
spin states. Indeed, in Fock space and for spin- 1/2 particles, the occupation of each site 
i can be characterized by the states In^j, > with n it0l G {0,1}, that is states 1 0, >, 
1 1, >, |0, 1 > and |1, 1 >. But in the case of single occupancy the states |0, > and |1, 1 > 
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which are excluded as unphysical in the present case have to be eliminated. This is done 

ft 

by means of a projection procedure proposed by Popov and Fedotov p and generalized to 



SU(N) symmetry in ref. 

Introducing the projection operator P = e l ^ N where N = J^fLfiv * s the number operator 
the partition function Z reads 



Z = Tr 



-pup 



where (3 is the inverse temperature. On each site i the contributions of states |0, > and 



11, 1 > to Z eliminate each other. Indeed 



< 0, 0| i e -/JH .e < 5* |0, >i + < 1, l^e^ 11 .e^* 2 \l, 1 >; 
+ < l,0\ l e~ l3H .e^\l,0> i + < 0,11^-^.6^10,1 >, t 
= i{< 1, 0| ie -^|l, >i + < 1, 0|<e-^|l, >i) 

Hence the partition function 



Z = -.Tr 

i 



(3) 



with the imaginary " chemical potential" /i 



2/3 



describes a system with strictly one particle 



per lattice site, in contrast with the usual method which introduces an average projection 
by means of a real Lagrange multiplier 



III. MEAN FIELD AND ONE-LOOP APPROXIMATIONS 

Following the usual procedure we transform the Heisenberg Hamiltonian into a bilinear 
fermionic expression using a Hubbard-Stratonovich decoupling. Starting from (3) this leads 
to 



Z 



Z n 



(4) 



where r is an imaginary time and 
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J i 

So[0(r)} = l^tfmr) -Bj) 

hi 

where (p stands for the Hubbard-Stratonovich decoupling fields and £ for the Grassmann 
variables. 

After integration over the bilinear fermionic terms which appear in the action Z 

takes the form 



Z = — I V0e~°' ' ; ! 



where 



S, 



eff 



/ drS [0{r)) - VZn2cfc£||&(u; = 0)|| +Tr{V-(G M 1 )"} 
Jo * 1 n 

v % n=l 



(5) 



and = 0) is the Fourier transform of 0i(r). The propagator Go and Mi are defined 

in matrix form in appendix IVI Al 

In a loop expansion beyond the mean- field approximation (p the effective action given by 
(0) is expanded in a Taylor series 



S, 



eff 



s. 



eff 



dS, 

+ 

[0] 80 



eff 



[0} r 2 dip 2 [0} 



S0 2 + O{50 3 



to second order (one-loop contribution) in the fluctuations Sip of = 0+50. Since the mean 

50 = one gets the set of coupled self-consistent 



field (p is chosen in such a way that 9S gl f 



equations 



E^ 1 



ipj 



B, 



2 0i 



which fixes the fields (p. 

In the following we consider a Neel mean-field order (pi(r) = (—l) 7T ' Ti z e z = 0\e z where 
7? is the Brioullin spin sublattice vector. The magnetic field applied to the system is also 
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chosen to be aligned along the direction e z . The partition function can be decomposed into 
a product of three terms 

Z = Zmf-Z zz .Z + - 
where Zmf, Z zz and are given by 



with 



Zmf 

Z zz 

z+_ 



~ S eff 

e : 



^zz 
^0 



Vip z e 



1 ° °eff 

2 d v z2 



5<p z 



z, 



1 ° °eff 



■I- / i _n 2 dip+d<p 

"VP > <P > e 



Sip+ .S<p 



s, 



eff 







J2ln2ch^\m 



'J 



1 d*S eff 



2 d^ z2 

1 d*S eff 



2 dtp+dtp- 











-i 



UJ l,J 



I) 



UJ l,J 



J; 



2 lJ 
2 J y 



(§# ) ) = 0) 



2 - vj3 J 13 



(6) 



Zmf is the mean field contribution, Z zz and Z + _ are the one-loop contributions respectively 
for the longitudinal part S(p z and the transverse parts of (p, S(f + ~ , which take account of the 
fluctuations around the mean- field value (p z . 

The contributions Z zz and Z + _ are quadratic in the field variables 5tp z , 5(p + ~ and can be 
worked out in the presence of a staggered magnetic field B z . Studies involving a uniform 
magnetic field acting on antiferromagnet quantum spin systems can also be found in ref . |s| . 
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IV. MAGNETIZATION AND SUSCEPTIBILITY OF D-DIMENSIONAL 

SYSTEMS 

A. Magnetization 



The fields {<fi} can be related to the magnetizations {mj as shown in appendix IVI Bl and 
the free energy can be expressed in terms of this order parameter, see appendix IVI CI The 



magnetization m on site i is the sum of a mean field contribution m = — -| 9ln ^ F , 
contribution <5m+_ = — \ dln ^,~ and a longitudinal contribution Sm zz — 1 : 



small magnetic field B a linear approximation leads to m = rh + <5m 2Z + <5m + _ where 



a transverse 
For a 



rn 



-th—D\J\m 
2 2 1 1 



<5m . 



E 



8mAm (1 — 4m 2 



keSBZ 



[l + 2D\J\Am ) 
4m 



- E 



4m 2 ) 2 



(l + 2D|J|Am (l 



1- 



keSBZ 



N p is the number of spin-1/2 sites, Am : 
appendix IVI CI for details of the derivation. 
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l-4.fn a 



th /3L>|J|m A /l-7 



l-f D|J|(l-4.m 2 ) 



and 7fc = X] cos(k.ff), see 



ridn.n. 



At low temperature (T — > 0) the magnetization goes over to the corresponding spin-wave 



expression 
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12L Il3l | , which reads 



m 



1 



- E 



1 



1 



I3D\J\ 
2 



l- 7 j 



1-7 2 



feessz — \ 2 \/ - < fc ; y - < fe 
Figure ^ shows the magnetization m in the mean-field, the one-loop and the spin wave 
approach for temperatures T < T c where T c = D\J\/2 corresponds to the critical point. 
One observes a sizable contribution of the quantum fluctuations generated by the loop 
contribution over the whole range of temperatures as well as an excellent and expected 
agreement between the quantum corrected and the spin wave result at very low temperatures. 

The magnetization shows a singularity in the neighbourhood of the critical point. This 
behaviour can be read from the analytical expressions of 5m + _ and Sm zz and is generated by 
the \k\ — mode which leads to jr — 1 and by cancellation of m. The Neel state mean- field 
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FIG. 1: Magnetization in a 3D Heisenberg antiferromagnet cubic lattice. Dashed line : Mean field 
magnetization, Dotted line : Spin wave magnetization, Full line : One-loop corrected magnetiza- 
tion. 

approximation is a realistic description at very low T. With increasing temperature this 
is no longer the case. The chosen ansatz breaks a symmetry whose effect is amplified as 
the temperature increases and leads to the well-known divergence disease observed close to 
T c . Hence if higher order contributions in the loop expansion cannot cure the singularity 
the Neel state antiferromagnetic ansatz does not describe the physical symmetries of the 
system at the mean- field level at temperatures in the neighbourhood of the critical point. 
Consequently it is not a pertinent mean-field approximation for the description of the system. 

The discrepancy can be quantified by means of the quantity where Am = m — fn = 
5m zz + 5m + _. Figure 2 shows the result. The relation < 1 (Ginzburg criterion) fixes a 
limit temperature Tu m above which the quantum fluctuations generate larger contributions 
than the mean-field. For 3D systems this leads to Tu m ~ 0.8T C , see figure EI 

The pathology is the stronger the smaller the space dimensionality. It is also easy to see 
on the expression of the magnetization that, as expected, the contributions of the quantum 
fluctuations decrease with increasing D. As can be seen in the figure the saddle point 
breaks down earlier in two than in three dimensions. 

In fact, the Heisenberg model spin wave spectrum shows a Goldstone mode as a conse- 
quence of the symmetry breaking by the Neel state. When |fc| goes to zero 
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FIG. 2: Ginzburg criterion -4^ for the 3D Heisenberg model. 



Am 

m 




FIG. 3: Comparison of the Ginzburg criterion for 2D (dashed line) and 3D (full line) Heisenberg 
model. | As | > 1 for T > at 2D and T > 0.8T C at 3D. 



Wjf = ZDSyJl—% 

lim u^ ~ |fc| 

fc^O 

The zero mode destroys the long range order in ID and 2D as expected from the Mermin- 
Wagner theorem jl^ . 
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In the case of the XXZ-model the Hamiltonian of the system can be written 
H XXZ = J2 ( S ? S * + S i S ] + (1 + 5 ) S * S i) 

<ij> 

where 5 governs the anisotropy. In this case the excitation spectrum shows a finite \k\ = 
energy lu^ 



j ^ 2 



lim luj: ~ 



1 




In appendix IVI Dl we develop explicitly the expressions of the free energy, magnetization 
and susceptibility By examination the expressions show that the zero momentum mode is 
no longer responsible for a breakdown of the saddle point procedure near T^ xz = D ^ J 2 +A ^ 
However the magnetization of the XXZ-model remains infinite near T* xz . This is due to 
the common disease shared with the Heisenberg model that the mean field magnetization 
appearing in the denominator of 6m + _ goes to zero near the critical temperature. One 
concludes that the mean-field Neel state solution makes only sense at low temperatures, 
that is for T < Tu m , whatever the degree of symmetry breaking induced by the mean-field 
ansatz. 



B. Susceptibility 

We consider the parallel susceptibility x\\ which characterizes a magnetic system on which 
a magnetic field is applied in the Oz direction. The expression of \\\ decomposes again into 
three contributions 



XII 



1 d 2 T 



XMF + Xzz + ^ 



B=0 



with 
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FIG. 4: Parallel magnetic susceptibility at 3D for Heisenberg model. Dotted line : Spin wave 
susceptibility. Dashed line : mean field susceptibility x\\ MF - Full line : total susceptibility (x\\ MF + 



MF 



-Am 



f (1 - 4m 2 ) 



X\\zz 



1 + §£>|J|(1 



/?g|J|7g 
2 



Af„/3 



E 

keSBZ 



Am 2 ) 

2 

Am 2 (1 + 4m 2 ) 



XII 



L E 



1 - 

1/3(1 



2 



(1 - 4m 2 ) 2 



2D|J|Am 



+ 



keSBZ 



2 
1 



sh 2 J\fh) 



sh 



^iJimyr^ 2 



o sh2/?L>|J|m 
| (1 - 2D| J|Am x0 ) 2 - /3 (D| J|Am x07A :) 2 



7e 



(3D\J\m^ I -7 



The behaviour of xy is shown in Figure 0] in which we compare the mean-field, spin wave 
and the one-loop corrected contributions for a system on a 3 D cubic lattice. One observes 
again a good agreement between the quantum corrected and the spin wave expressions at 
low temperatures. For higher temperatures the curves depart as expected. The mean-field 
contribution remains in qualitative agreement with the total contribution. 
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V. CONCLUSION 

In the present work we showed the contribution of quantum fluctuations in the description 
of quantum spin systems at finite temperature governed by Heisenberg Hamiltonians in space 
dimension D. The mean-field ansatz is taken as a Neel state. The number of particles per 
site is fixed by means of a rigorous constraint implemented in the partition function. It has 



been shown elsewhere 15[ that this fact introduces a large shift of the critical temperature 
compared to the case where the constraint in generated through an ordinary Lagrange 
multiplier term. 

At low temperature the magnetization and the magnetic susceptibility are close to the spin 
wave value as expected, also in agreement with former work 7}. The quantum corrections 
are sizable even for low temperatures. They increase with increasing temperature. 

At higher temperature the quantum contributions grow to a singularity in the neighbour- 
hood of the critical temperature. The assumption that the Neel mean-field contributes for 
a major part to the magnetization and the susceptibility is no longer valid. Approaching 
T c the mean field contribution to the magnetization goes to zero and strong diverging fluc- 
tuations are generated at the one-loop order. This behaviour is common to the Heisenberg 
and XXZ magnetization. In addition the Neel order breaks SU(2) symmetry of the Heisen- 
berg Hamiltonian inducing low momentum fluctuations near T c which is not the case in the 
XXZ-model. 

The influence of quantum fluctuations decreases with the dimension D of the system due 
to the expected fact that the mean-field contribution increases relatively to the one-loop 
contribution. 

For dimension D = 2 the magnetization verifies the Mermin and Wagner theorem ^| 
for T^O, the fluctuations are larger than the mean field contribution for any temperature. 
Therefore another mean-field ansatz has to be found in order to describe the correct physics 
at finite and not too low temperature. 



The authors would like to thank Drs. D.Cabra and T. Vekua for instructive discussions. 



12 



VI. APPENDIX 



A. Matrices G and Mi 



After integration over the fermionic degrees of freedom in equation (J1J) the partition 
function takes the form 



Z = — f V(pe~ { f° ' *rSo[v(T)]-lndeHpM\} 

Zo J 

z Q J * 



where 



Mi 



ipS P , q + \^\{p ~ q) \<Pi (p - q) 

\vt (p - - i^(p - ?) 

p and 9 refer to modified fermionic Matsubara frequencies, p — Uf — \i — ^(n+ 1/4) and n 
is an integer, see jf|. M can be put in the form 



where 
G = 



M 



det 



fez Mm - iv? (p - e = °)<Wl 



1 1,^ 



det G„ 2 



1 1 



<pt(p-q = o)5, 



P>9 



ciet 



Mi 



5^(P-?) 5^t(p-?) 

-l s f!(p-q) 



with 50i(p — q) — 0i{p — q) — <Pi(p — q — 0)5 P , q - The second term in the expression of M 
corresponds to the quantum contributions. The expression In det (/3M) can be developed 
into a series 



lndet{(3M) = In det (3 [-G \l - G M ly 

= Indeti-pGo 1 ) +Trln(l - GqM!) 

00 ^ 

= Indeti-PGo 1 ) -Tr{^-(G Mi) n } 



n=l 
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The first term lndet(—(3G 1 ) leads to the expression £V ln2ch^\\<Pi(u = 0)||. The first term 
in the sum gives the contributions at the one-loop level. 

B. Relation between the Hubbard-Stratonovich mean fields (pi and the mean-field 

magnetizations m,j. 

Using rhi = — 9 q^ z f one gets 

2 ,_i 
(fj = —th Ivrii 



P 



2 

-th~ 1 2m i = + Jij^j 



P 

j 

Starting from the Neel state rhi = (— l)' l fn + (-l)'Am(B). For a weak magnetic field 
the second term can be linearized Arh(B) = Am B. The coupling strength matrix acting 
between nearest neighbour sites is taken as .1^ = J Yl <Kn ~~ =t V) with J < 0. Then 

ffen.n. 

rhi = (-iym + (-iyAm(B) 

= [(-iyB + Z\J\(-iy(m + Am(B))} 

(-l) i (m + Am{B)) = \th^{-\y [B + 2.D.\ J\.(m + Am(B))} 
By means of a Taylor expansion around B = : 



(-iy(m + Am(B)) = -th-(-iy.2.D\J\m 

+ l^(-iy[l + 2.D\J\Arh ] 



l-th 2 ^{-l) l 2D\J\m 



.B 



+ 0{B 2 ) 

By identification one gets Am — § [1 + 2.D| J|Am ] [1 — 4m 2 ] and finally 
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Am 



m = -th—D\J\m 
2 2 1 1 

Arh(B) = Arh .B 
f(l- 4m 2 ) 



1 - §D\J\ (1 -4m 2 ) 
where D is the lattice dimension. 



C. The free energy and the terms 6cp z and 5ip + in equation 1^1 



ld 2 S, 



<ff 



Substituting m; = (— l) l (m + Am(B) in ^ d zi 
tion |H1 leads to 



dip ana 2dv+dip - 



5ip + 5ip of equa- 



& 



(l-4(m + Am(B)) 2 ) 



(2) 

[2aL 
[26]. 



[M u + (-i) c 

(m + Am(B)).(B + 2D\J\(m + Am(B))) 
\{B + 2D\J\(m + Am(B))) 2 + cj 2 ] 
iu;(m + Arh(B)) 
[(B + 2D\J\(m + Am(B))) 2 + u 2 } 



Integrating out the fluctuations dp away from the mean field p z leads to 
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Am 



m = -th—D\J\m 

Am(B) = Am .B 
§ (1-A.m 2 ) 



1 - §£>|J| (1-A.m 2 ) 
Tmf =N p D\J\ (m + Am(B)f - ^ln ch [B + 2D\J\(m + Am(B))]j 

S ^=^ E in[l-fffl^) 2 [l-4(m + Am(S)) 2 ] : 



keSBZ 



/Smh(l[[B + 2D\J\(m + Am{B))] 2 - [2D\ J\ lu {m + Arh(B))\ 

= 4 Y In ^ 775 ; 

P f-i V Sinh(f[B + 2D|J|(m + Am(S))] 



1/2 



D. The XXZ-model 



The XXZ Hamiltonian 



H XXZ = ~ J - J2 ( S i S ! + S i S J + (! + 



2 

<ij> 



leads to a critical temperature 



xxz D\J + A\ 



2 

A = J 5 



and a mean magnetization 



m = hh^D\J + A\m 



£(1-4 m 2 ) 

Am(B) = Am .B = B.- ' - 



1 - %D\J + A| (1 -4.m 2 ) 
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1. Free energy 



J~mf =N p D\J + A\ (m + Am(B)) 2 - ^ln ch [B + 2D| J + A|(m + Am(B))]j 
8T ZZ =i J! '» " ( ^ |J 2 +A ' 7g ) 2 [1 - 4(m + Am(D)) 2 ] : 



keSBZ 



*r+-=± E M 



/Sinh f § ([5 + 2D | J + A|(m + Am(D))] 2 - [2D|J| 7jE (m + Am(B))] 2 ) 



1/2 



keSBZ 



Sinh (f [5 + 2D | J + A|(m + Am(D))]) 



2. Magnetization 



m = 



m = 



Sm,r. = 



m + 5m zz + 5m+- 
1 dTuF 



N p OB 
1 dT zz 



N p dB 



B=0 



B=0 



N p (3 



E 

keSBZ 



8.mAm (l-4m 2 ) (W±^y 



1 - 



//3£>|J+A| 7 g 



5m + _ = 



\ 2 
(1 + 2D|J + A|Am ) 



) 2 (1 -4m 2 ) 



4m 



- E 

AT, ^ 



1 + 2DIJ + A|Am (l - 



J+A 'fc, 



keSBZ 



\l l -{lh^Y th (pDlJlfhy/l-^-yz)' 



3. Susceptibility 
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X\\ 



MF 



XL 



X\\. 



1 d 2 T 



N p dB 2 



Xmf + Xzz + Xh 



B=0 



-Am 



xo 



§(l-4.m 2 ) 
1 + §D\J + A| (1 -4m 2 ) 



N p f3 ^ 



^4^) Am 2 (l + 4m 2 ) 



keSBZ 



/3D|J+A| 7 r 



) (l-4m 2 ) 2 



- E 



1/3(1 -2DI J + A| Am 



a; 
+ 



, 2 s h 2 (/?£>[ J + Aim) 



sh 2 (^|J + A|m A /l-(^) 2 ) . 

fl / n irlA , 2 ^2/?D|J + A|m^l-(^ 7g ) 

- /5 (£>K I Am x o7fe) / 

/3 J D|J + A|m A /l-(^) 



^(l-2/J|J + A|Am x0 ) 2 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8] 

[9] 
[10] 

[11] 
[12] 

[13] 



T. Senthil, A. Vishwanath, L. Balents, S. Sachdev, M. Fisher, Science 303 1490 (2004) 
T. Senthil and M. P. A. Fisher, |cond-mat/0510459 



P. Ghaemi, T. Senthil, lcond^mat/0509073 (2005) 



T. Morinari, cond-mat/0508251 (2005) 

P.A. Lee, N. Nagaosa, X.-G. Wen, |cond-mat/0410445] (2004) 

V. N. Popov and S. A. Fedotov, Sov. Phys. JETP 67 535 (1988) 

S. Azakov, M. Dilaver, A.M. Otzas, Eur. Phys. J. B, 22 53 (2001) 

M. Kiselev, H. Feldmann, R. Oppermann, Eur. Phys. J. B 22 53 (2001) 

A. Auerbach, Interacting electrons and quantum magnetism, Springer- Verlag, 1994. 

A. Auerbach and A. Arovas, Phys. Rev. B38 316 (1988) 

J-I.Igarashi, Phys.Rev. B, 46 10763 (1992) 

E.Manousakis, Rev.Mod.Phys., 63 (1991) 

T. Holstein, H. Primakoff, Phys.Rev., 58 1098 (1940) 



18 



[14] N.D. Mermin, H. Wagner, Phys.Rev.Lett. 17, 1133 (1966) 
[15] R. Dillenschneider and J. Richert, cond-m at / 0509724 



